Response to Greiter's Comment 
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Greiter |l| claimed erroneously that the n excitation |^| of the Hubbard model has an energy of 
the order of U. This mistake originates from a inconsistent treatment of the two particle Hartree 
energy and the Hartee correction to the chemical potential. In any self-consistent calculations, these 
two contributions cancel exactly, as shown by Kohn, Luttinger and Ward We also show that his 
interpretation of the finite-size studies is inconsistent. 
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In ref , we studied the n resonance within the Hub- 
bard and the t — J model. We calculated the spectrum 
of the 7r operator 
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within the T matrix approximation (which is equiv- 
alent to the Hartree-Fock approximation in the two 
particle channel), and showed that within this approx- 
imation, the spectrum consists of a sharp pole at 
uiq = — 2/i + J/2 (1 — n), where \x is the chemical poten- 
tial approperiate for the t-J model. Meixner et al. Q 
have checked this result within exact numerical diago- 
nalization, and found excellent agreement. 

Greiter jl| claimed to have found an error in our 
Hartree-Fock calculation, and claimed to have shown 
that the energy spectrum of the 7r operator is of the or- 
der of U, and it is incoherent. Both of these claims are 
wrong, and we shall trace his mistake in this note. 

Before proceeding to the calculations, let us first show 
the physical origin of his error. The n operator inserts 
two extra particles in triplet channel into a fermi system. 
In the grand canonical ensemble, the energy of the n op- 
erator is defined with respect to twice the chemical po- 
tential. Since the ir operator is a spin triplet, the on-site 
Hubbard U interaction can not enter the multiple scat- 
tering process between these two extra particles directly 
but only through a Heisenberg exchange term J. Within 
the Cooper approximation, the only role of the fermi sea 
is to block the filled states in the multiple scattering pro- 
cess, and the chemical potential is simply given by the 
fermi energy of the free fermi sea. At this zeroth order, 
one can clearly see that the energy of the tt operator is 
independent of U. To the first order in U , one can con- 
sider the interaction between the inserted electron with 
the rest of the electrons inside the fermi sea, and there is 
a non-zero Hartree contribution 2£/?ij. However, to this 
order, there is also a Hartree contribution to the chemical 
potential. Since the energy of inserting a single electron 
into the fermi sea is also raised by Un^, these two con- 
tributions cancel each other exactly, and therefore the 



energy of the 7r operator is unchanged with the inclusion 
of the Hartree process. There is no Fock process for the 
Hubbard U. Greiter reached the erroneous conclusion 
because he only took the Hartree contribution to the two 
particle energy, while overlooking the correction to the 
chemical potential which is of the same order. 

In a series of classic papers, Kohn and Luttinger, Lut- 
tinger and Ward || have shown that this type of can- 
celation occurs order by order in perturbation theory. 
Therefore, one can make a stronger argument, namely 
the energy of the 7r resonance will remain finite in the 
large U limit, if the perturbation series converges in this 
limit. However, this assumption may break down, and 
the Hartree treatment of the interaction would certainly 
be inadequate. Unfortunately, no systematic calculations 
of the 7T spectrum have been carried out in the limit where 
double occupancy is strictly forbidden. Greiter's com- 
ment is purely based on the Hartree-Fock approximation, 
and did not add anything to the analytical calculations 
in this limit. Meixner et al [Q presented strong numerical 
evidence that the tt spectrum remains sharp in this limit, 
and have shown furthermore the low energy of the peak 
scales inversely with U. Greiter tried to argue that the 
7T peak could be due to mixing with the spin correlation 
function. Close to half-filling, there is considerable pair- 
ing fluctuation, and we indeed first suggested that these 
two channels are mixed. If one only looks at the corre- 
lation functions close to half filling, there is an obvious 
"chicken and egg" problem in interpretation. However, 
Meixner et al. found that the sharpness of the 7r peak 
persists over a wide doping range, up to n = 0.4 0. In 
this doping regime, there is no sign of any pairing fluctu- 
ations, and the sharpness of the tt peak must be intrinsic, 
completely decoupled from the spin correlation function. 

For most readers familiar with the basics of the 
Hartree-Fock or T matrix approximation, the above dis- 
cussion should have adequately addressed the elementary 
nature of the debate raised by Greiter and the source 
of his error. However, in order to make this note self- 
contained, we shall show the details of the calculations 
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in the remainder of this note. 
The Hubbard model is given by 



U J2 



(2) 
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For a given density of electrons n, the chemical potential 
/i is implicitly determined from n = —dEo/dfi, where Eq 
is the ground state energy. Order by order in pertur- 
bation theory, the dependence of Eq on /i changes. In 
order to keep the density of the system fixed at any or- 
der, one has to renormalize [i consistently. In general, fj, 
determined from the above relation changes from order 
to order in perturbation theory. This point is discussed 
in detail by Kohn, Luttinger and Ward ||. 

Let us see how the Hartree-Fock approximation renor- 
malizes the chemical potential. For non-interaction elec- 
trons n — ^ fc 9(fj, — 6k)- This determines the chemical 
potential in zeroth order to be the non-interacting fermi 
energy /l*o = cf- Next we can consider the commutator 
between the Hamiltonian and the one electron creation 
operator, evaluated within the Hartree-Fock approxima- 
tion, 
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To this order, the ground state wave function is modified 
to be 



cLl°> = ° »/ ( k - M + Un-„ < 
Cfc<x|^) = if e k - fj, + Un^ a > 
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In this new ground state, the density of electrons is given 
by 

11 = 4 £< n l c L<*r|n) = 4 £ 6 0* - Un x - e k ) (5) 
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Comparing this expression with the density of the non- 
interacting electron gas, we see that in order to fix the 
density, the chemical potential at this order is given by 
ji = Uni + ep. 

At the same order in U, the energy of the tt operator is 
determined by evaluating the commutator in the Hartree- 
Fock approximation, 

[H, Tr' d ] = 2(Uni — ^)tt\ + (terms of the order of J) (6) 

Here the terms of the order of J arises from multiple 
scattering process of the added particles with each other. 
In lowest order, the J terms can only enter if we add a 
Heisenberg exchange term to the Hubbard model. This 
term is introduced in order to be able to show a simple 
explicit result within the lowest order Hartree-Fock the- 
ory. Its introduction is not necessary if one performs a 
full diagrammatic calculation of the Hubbard model, as 
we will show later in this note. Anyway, the effect of J is 



not challenged by Greiter's comment and we will not dis- 
cuss it explicitly. Since the added electrons form a spin 
triplet, the on-site U does not contribute to their mutual 
interaction. The first term arises from the Hartree in- 
teraction of the added electrons with the electrons inside 
the fermi sea and represents a self-energy correction to 
a one electron's Green's function. At this point, Greiter 
argued incorrectly that the energy of the tt operator is 
of the order of U. Because of the renormalization of the 
chemical potential discussed in the previous paragraph, 
we see that the U term cancel and 2(Z7nj — fi) simply 
reduces to — 2ej?, which is manifestly independent of U. 

The same effect can of course also be addressed di- 
rectly in the canonical formalism, without introducing 
the chemical potential. In this approach, the energy of 
the 7r operator is given by 

c^o = [E W (N + 2) - Eo(N)] - 2[E Q (N + 1) - E (N)} (7) 

where E 7T (N + 2) is the energy of a lowest excited state 
of the N + 2 particle system with total spin one and 
momentum Q. The numerical calculations of Meixner et 
al. H| uses this definition for the energy of the tt operator. 
If we neglect the interaction between the added electrons 
with the background, the first term is of the order of j[] 
and the second term is simply —2ep. If we include the 
Hartree interaction with the background, both terms are 
raised by 2J7n;, and cancel each other. We obtain the 
same expression as before. 

The reader may wonder if this type of cancelation hap- 
pens only accidentally in lowest order of perturbation 
theory. In fact, by using a powerful theorem due to Lut- 
tinger and Ward |[| , one can show that it occurs at every 
order in perturbation theory. The energy of the 7r reso- 
nance within the pure Hubbard model can be obtained 
by solving the following exact integral equation 



(8) 



Here T lrr is the sum of all two particle graphs which do 
not fall into two separate pieces when one cuts two Greens 
function lines, and G is the exact Greens function given 
by 
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e p - fi- 



(9) 



* For the tight binding Hamiltonians the dispersion is given 
by e p = —2t(cosp x + cosp y ), and adding two particles with a 
center of mass momentum Q — (tv,tt) results in a zero kinetic 
energy e p + e- P +q = 0. The first term in ([?]) is therefore 
determined by a repulsive interaction of two particles in a 
triplet state at nearest sites 
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and £(w,p) is the exact self-energy. Because of the spin 
triplet nature of the n resonance, the symmetry of the 
vertex function exactly projects out the on-site interac- 
tion part in T lrr . The interaction on the nearest neighbor 
sites due to the remaining part in T lrr can be defined to 
be the effective spin exchange constant J. £(w,j?) can 
involve many corrections, including the on-site interac- 
tion, just as we have seen in the first order calculation. 
However, the chemical potential also receives a similar 
correction. Luttinger and Ward (^] obtained an explicit 
formula relating these corrections, 

e F = fx - T,((x,p F ) (10) 

which is valid to all orders of perturbation theory. We 
can identify the constant part of £(w,f>) with the correc- 
tions due to an effective instantaneous on-site interaction. 
From the Luttinger Ward formula, we see that any such 
correction is exactly canceled by a corresponding shift in 
the chemical potential. The remaining part of the self- 
energy renormalizes the bare electron into a dressed elec- 
tron with modified fermi parameters. It is plausible that 
in the low energy sector, the above integral equation has 
a form similar to the lowest order Hartree-Fock equation 
and a low energy collective pole exists as its solution. 

In conclusion we have identified the error in Greiter's 
comment: it results from a inconsistent treatment of 
the Hartree correction to the two particle energy and 
the corresponding correction to the chemical potential. 
We demonstrated explicitly the cancelation of these two 
types of corrections within the lowest order Hartree-Fock 
or the T matrix calculation. Using a theorem due to Lut- 
tinger and Ward, we argued that this type of cancella- 
tion occurs to all orders of perturbation theory. Greiter's 
interpretation of the numerical results overlook a large 
data set on the 7r resonance. In the parameter regime 
where pairing correlations are present, we argued previ- 
ously that the spin correlation and the ir correlation may 
mix. However, tt resonance has been clearly identified nu- 
merically in the regime of low electron density where the 
pairing correlation is absent, and is therefore intrinsic. 
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